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Abstract. We examine the possibility that pertinent impurities in a condensed 
matter system may help in designing quantum devices with enhanced coherent 
behaviors. For this purpose, we analyze a field theory model describing Y- shaped 
superconducting Josephson networks. We show that a new finite coupling stable 
infrared fixed point emerges in its phase diagram; we then explicitly evidence that, 
when engineered to operate near by this new fixed point, Y-shaped networks support 
two-level quantum systems, for which the entanglement with the environment is 
frustrated. We briefiy address the potential relevance of this result for engineering 
finite-size superconducting devices with enhanced quantum coherence. Our approach 
uses boundary conformal field theory since it naturally allows for a field-theoretical 
treatment of the phase slips (instantons), describing the quantum tunneling between 
degenerate levels. 
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Figure 1. Sketch of the YJJN; Inset: graphical exact solutions for the energy levels 
at 5 = 9/8. 

For engineering quantum devices one has often to tame the decoherence arising 
from the interaction of a pertinent two-level system with both the control circuitry and 
the quantum modes lying outside the subspace spanned by the two operating states. 
An important source of decoherence arises when the total state of the two-level system 
and of its environment evolve towards an entangled state. If a system is coupled to 
more than one bath, and its entanglement with each one of the baths is suppressed 
by the other (s), decoherence may be frustrated [HE]- In this paper, we evidence how 
frustration of decoherence may arise from the existence of a finite coupling fixed point 
(FFP) in the phase diagram of the quantum theory describing the device. 

Existence of finite coupling fixed points in condensed matter is a rare instance 
realized, so far, only in quantum systems with pertinent impurities. Remarkable 
examples of systems exhibiting attractive FFP's are provided by the two-channel single- 
impurity [3] and two-impurity [1] overscreened Kondo models, as well as by F-shaped 
junction of quantum wires [5]. At variance, F-shaped junctions of one- dimensional 
atomic condenstates [6] exhibit a repulsive FFP, signaling the existence of a new 
transition point between stable weakly and the strongly coupled phases. 

Boundary conformal field theories [7j are a natural setting to investigate stable 
phases and phase transitions of quantum impurity systems, once the quantum impurity 
is traded [7| for a boundary interaction, involving only a subset of the bulk degrees of 
freedom: the boundary interaction is then renormalized by the bulk degrees of freedom, 
and the infrared (IR) behavior is determined by the stable fixed point (s) in the phase 
diagram. 

Superconducting Josephson devices are not only promising candidates for reahzing 
quantum coherent two-level systems [8j, but also provide remarkable realizations of 
quantum systems with impurities, whose phase diagrams, in the simple cases so far 
investigated, admit only two fixed points: an unstable weak coupling fixed point (WFP), 
and a stable one at strong coupling (SFP) [9]. The approach developed in Ref.[9] 
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naturally allows for a field-theoretical treatment of the phase slips describing quantum 
tunneling between degenerate levels, and provides remarkable analogies to models of 
quantum Brownian motion on frustrated planar lattices fiUl [TT] . When an effective two- 
level quantum system is operated near by the WFP or the SFP, there is no frustration 
of decoherence, since, at strong coupling, there is not even quantum tunneling between 
the degenerate states while, at weak coupling, there is full entanglement between the 
two degenerate states and the plasmon modes. In the following, we shall show that a 
FFP emerges in a y-shaped Josephson junction network (YJJN), and that it may be 
pertinently used to engineer two-level systems with enhanced quantum coherence. 

A YJJN is realized by joining a circular Josephson junction array C to three finite 
Josephson chains via weak links of nominal strength A (see Fig{T]). C is pierced by a 
dimensionless magnetic flux and is joined to one of the endpoints of the three chains 
(inner boundary); the other endpoints (outer boundary) are connected to three bulk 
superconductors at fixed phases ipi {i = 1,2,3). For simplicity, we assume that all the 
junctions in the YJJN are of strength Ej and that A -C Ej. The Hamiltonian describing 
C is given by 



He 



where A0, 



d 



—%- 



2Ej E cos 

i=l 



(1) 



is the phase of the superconducting order parameter at 



grain and Wg is a gate voltage, li Ej / E^. <^1^ W^^'^ = + | + /i, with integer and 
< /i < 1/2, the low-energy dynamics is governed only by the two states with total 
charge equal to A^ and to A^ + 1. 

The procedure outlined in Ref.[9] allows to describe the three finite chains with a 
Tomonaga-Luttinger Hamiltonian 



-E 



dx 



1 (d^,. 



dt 



dx 



(2) 



In Eq.(l2]) describe the plasmon modes of the chains, and g and v depend on the 
constructive parameters of the network [H]. 

Fixing the phase at the outer boundary of the chains sets Dirichlet boundary 
conditions on $i(x) aX x = L: ^i{L) = (fi. Since we require that the charge tunneling 
between C and the inner boundary of the three chains is described by a Josephson- 
like interaction, with nominal strength A -C Ej, one should use Neumann boundary 
conditions at the inner boundary, i.e. = Vi. This allows to write the tunneling 

Hamiltonian as Ht = — AX)f=i cos[$i(0) — (po^]- 

A boundary field theory approach allows to trade He + Ht with an effective 
boundary Hamiltonian, Hb, involving only $i(0), and given by 



Hb = -2Ew E • cos[«. ■ x(0) + 7] : 



(3) 



i=l 
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with = -^[^l{x) - ^2{x)l X2{X) = + $2(x) - 2<l>3(x)], «i = (1,0), 

1_ 

«2 = (-i^), «3 = (-i-^), 7 = tan-i[3tan(f)], and Ew = (f ) ' ^vk, with 



24(£; )2^2 y 1 + 2sin (-3). The colons : : denote normal ordering with respect to 
the ground state of the plasmon modes, |{0}). In the following, we shall argue that, 
for 7 = 7r/3, there is a finite range of values of g, for which a YJJN supports a FFP: 
this results from the fact that, for this value of 7, the two plasmon baths Xi ^''^^ X2, 
cooperate to destabilize both the SFP and the WFP. 

The perturbative second-order renormalization group (RG) equation for the running 
coupling strength G = LEw, given by 



dG 




1--]G-2G^ , (4) 



shows that Hf, is a relevant perturbation for g > 1, while it is irrelevant for g < 1. In 
Eq.(jl]), Lq is a pertinent reference length scale. The strongly coupled fixed point (SFP) 
is reached when the running coupling constant G goes to 00. The fields Xj{x), j = 1, 2, 
now obey Dirichlet boundary conditions at x = and Xi(0), X2(0) are determined by the 
manifold of the minima of the effective boundary potential (Eq.Q). One sees that for 
< 7 < 7r/3, the minima lie on the triangular sublattice A, defined by (xi(0), X2(0)) = 
(27rmi2, ■^[2'7rmi3 + '7rmi2]), while, for n/S < 7 < 27r/3, the minima lie on the triangular 
sublattice B, given by (xi(0), X2(0)) = (27rmi2-^, ;^[27rmi3+7rmi2-vr]), with mi2, mia 
relative integers. From Eq.([3]), one sees also that the difference in energy between 
the sets of the minima forming the A and B sublattices is given by ~ -E'jsin(7 — |). 
The manifold of the minima is depicted in Figj2l where the instanton connecting the 
degenerate minima of the honeycomb lattice emerging when the A and B sublattices 
are degenerate, is shown. 

Following the approach outlined in Ref.[lO], instanton effects may be taken into 
account through 

Hb = -Yj2{T-V,{0)+T^mO)} . (5) 

i=l 

In Eq.(l5]), T is an effective isospin operator, connecting two neighboring minima of the 



[+)f«i-e 



with 



honeycomb lattice of the zero-mode eigenvalues, Vi{Vi) =: exp 
being the dual fields of x, while Hb is the "dual" boundary Hamiltonian of Eq.Q. Y 
is an effective coupling defined as Y = Ej — Ew [12|. From the O.P.E. of the vertex 
operators entering Hb, the RG equation for the running coupling strength y = LY is 

dy A 4:g\ 2g 



For 7 = 7r/3 and 1 < (7 < |, neither the WFP, or the SFP, are stable. Accordingly, a 
minimal hypothesis for the phase diagram requires a FFP a.t y = y^,, with y^, finite. For 

instance, for (? = | — e, with e -C 1, one obtains y* ~ (^|^ ^ ^/e. 
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Figure 2. Some points on the two lattices A and B: at 7 = 7r/3, the two lattices 
are degenerate. In this case, the set of the minima of the boundary potential spans 
a honeycomb lattice, whose sites are connected by instanton trajectories as the one 
drawn in the figure. 



The energy of the minima may be varied by changing the phases of the three bulk 
superconductors since the eigenvalues (^1,^2) of the zero-modes of the fields Xj (it obeys 
Dirichlet b.c!) depend on the phases ipj as 



(Pl,P2) 



on sublattice A, and 



(7) 



B. In Eqs. fl71l8|l (^1.1,^2) are integers, Pi 



+ 2V2f32)) , (8) 
{(pi - Lp2)/{27ry/2), and 



on sublattice 

(32 = {ipi + ip2-2ip3)/{27rVQ). 

As it happens with other superconducting systems [8], also a YJJN supports a 
two-level quantum system, operating between two pertinently selected quantum states. 
Indeed, for ?/ -C 1 and near the SFP, the low-energy spectrum is given hj E = ^[p\'^+E', 
where p = {pi,P2) labels the zero-mode contribution, while E' comes from the plasmon 
modes: thus, for 7 = 7r/3, a pertinent tuning of Pi and P2 renders degenerate the zero- 
mode contributions to the total energy coming from two nearest neighboring sites of the 
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honeycomb lattice resulting from the degeneracy of the A and B sublattices (see FigJ2]). 
This happens, for instance, if Pi = 1/3^/2,^2 = 0: the two degenerate quantum states 
I I) and I I) -labelled by (^1,712) = (0,0) on sublattice A and by (ni,n2) = (1,0) on 
sublattice B- are macroscopically characterized by the opposite values of the Josephson 
current flowing across chain-1 and chain-2, namely: Ji = —I2 = ifc ^|^^ , I3 = 0. 

Quantum tunneling between the degenerate states is induced by Hb, with matrix 
element —Y. Setting P2 = 0, and Pi = l/3\/2 + 5/(2tx), with 5/2tx <^ 1, one easily gets 
an effective Hamiltonian for the two-level quantum system as 

H2 = tQ{5)\ + t{5)a' -Y(y+Vi{Q)-Ya-Vi{^) . (9) 

In Eq.(l9]) eo('^) = f (| + ^(^) ~ 3"^' ^'^'^ Pauli matrices, 5 is a control 

parameter determined by the phases and —Y[a^Vi{Q) + cr~Vi(0)] describes the 

interaction of the two- level system with the phase slip operators introduced in Eq.([5]) 

In the spin Hamiltonian describing the two level system in Eq.([n]), one sees that 
there is a z-component proportional to e((5), as well as an x-component proportional 
to Y . While e{5) does not get renormalized by the interaction with the two plasmon 
fields, Y is renormalized and its value measures the amount of entanglement between the 
two-level system and the plasmon modes bath. In particular, if Y is irrelevant, the two- 
level system decouples from the environment and behaves as a classical (Ising-like) spin, 
pointing along z. When this happens, no energy is dissipated into the environment, 
and the spectrum of the Hamiltonian in Eq.(l9]) is given by two classical states with 
oj = ±e(5). If y — s> 00, the effective field acting on the two-level system would again 
make it behave as a classical spin, pointing in the a;-direction: now, all the energy 
is dissipated into the environment and the spectrum of Eq.(l9]) is given by only an 
overdamped mode at = 0. Only when Y takes a finite value y^,, the competition 
between e{6) and may lead to the emergenge of the frustration of the decoherence 
of the two-level system, since there is the possibility that, for a pertinent choice of the 
control parameter e{6), there is a finite damping, resulting in two broad modes, centered 
at pertinent renormalized energies. 

To evidence the frustration of decoherence around the FFP, we compute the spectral 
density of the Hamiltonian in Eq.(l9]), given by x^(a;)/a;, where xli^) is the imaginary 
part of the transverse dynamical spin susceptibility p!]. The diagrams contributing to 
X±{.^) are shown in Fig.(l3]b): x±{^) is computed as a loop defined by the | |)-state 
propagating forward in (imaginary) time, and by the | |)-state propagating backward. 
It is given by 

/CO 
7ri9l^{-z)gu{^ + uj)+gl^{-z)g^^^{z + uj)} , (10) 

where gaai^) is the Fourier tranform of the propagator of the "spin" eigenstate \a) 

(^=T,i). 

For 7 = 7r/3 and g > 9/4, the boundary interaction is irrelevant, and one may 
neglect corrections to the amplitudes of order F^. This amounts to substituting gcrai^) 
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Figure 3. a): Graphical representation for the Schwinger-Dyson equations for 
g„„i{ijj); b): The "bubble" yielding the dynamical spin susceptibility; c): Graphical 
representation of the RPA summation implemented to compute X-l('^) near the FFP. 



with its noninteracting limit, g^\uj 



l/[iuj + cre((5)], yielding 
oc [5(cu + 2e(5)) - ^(cu - 2e((5))]/cu . 



'Ill 



Eq. lfTTl) shows that there is no entanglement (for g > 9/4) between the two- level 
quantum system and the plasmon modes. Since, in this range of g, there is no tunnel 
splitting between the two degenerate states, the system is classical and no quantum 
coherence emerges. 

For 7 = 7r/3 and g < 1, instantons provide a relevant perturbation and the IR 
behavior of the system is driven by the WFP. To compute x_l(co'), one now needs 
to substitute gjj\u}) with the dressed propagator, gf,.a-{to), drawn in Fig.([3]a), where 
the solid heavy line represents the fully dressed propagator, while the solid light line 
represents g^\uj). As a result, g„„i{uj) is given by 



(o)i-i, 



-\u) + Y^V,{u)}{[gf\ 



Y^T^{uj)} + Y^ 



with cr = — cr, and T(j{uj) is the Fourier transform of the propagation function 



({0}| : e=^^i®(^i) :: e'F^i®^^^^ 



|{0}) 



(12) 



^(l3) 



at frequency uo — ae{5). X_l(co') near the WFP is computed by taking the large-y limit 
of Eq.dlSl), yielding 



oc [|2e(5) + uj\^--^ - \lt[b) - 



(14) 
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Eg. (1141) shows that a large part of the spectral weight has moved now from the side 
peaks towards uj = 0, signaling the strong decoherence of the two-level system described 
by Eq.®. 

For 7 = 7r/3 and 1 < (7 < |, the IR behavior is driven by the FFP. A closed- 
form computation of x±{^) is now possible only for special values of g. For instance, if 
g = ^ — € with e -C 1, y^, is <^ 1, and one may compute x±{.^) by resorting to a RPA 
summation, graphically sketched in Fig.([3]c)). The result is 

1 

[X±]RPA(t^) 



io - 2y'e2(5) + Y^- Y^T[-1 - |7](-^) 
+ _ .L. . B,. • (15) 



uj + 2Je^{5) + r2 _ Y^T[-1 - h]{-uj 
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When writing [x±]rpa(i^) as a function of the dimensionless variable taking 
into account that the dimensionless variable y = YL — ^ as L ^ 00, one gets 

1 1 

[x±]rpa{x) oc 3 3 + o 3 , (16) 

X — 2A + e s^x^gT x + 2A -|- e a^x^s^ 

where A = ^^/WTY^. The imaginary part of Eq. f|T6|l shows two peaks centered 



around ±2y [e(5)]^ + (^y*)^, where y^, is the finite fixed point value of the running 
coupling constant. In FigJU we provide the plot of x "(co')/co'. 

A very special situation arises for 7 = 7r/3 when g = 9/8 since, for this value of g, 
the scaling dimension of the relevant instanton operators equals 1/2, just as it happens 



with fermionic operators. Indeed, for g = 9/8, the plasmon field : exp —■^iy/2gQ(x 



becomes a fermionic operator iIj{x + ivr) {—L < x < L), and the spin-1/2 operators 
may be fermionized according to 



a' — ^a^a-- , a+ — > a^e' , (17) 

where the zero-mode operator Pi ensures, for g = 9/8, the correct anticommutation 
relations between ip and the operators in Eq.( JT71) . As a result, the two-level Hamiltonian 
Eq.([9]) becomes 



i/per = -tv r dxV'^(x)^!^+e(5)(ata--)-rV^(0)ate^"V«'''_h.c.,(l^ 

J-L OX 2 

with twisted boundary condition 
ip{L) = exp 



H-L) , (19) 



where pi is the eigenvalue of the zero- mode operator Pi. A similar situation arises in 
the analysis of a spin-1/2 Kondo system at the Toulouse point [13]. In particular, the 
Hamiltonian in Eq. (fT8|l has been recently proposed [13] to describe two qubits at the 
end of a finite length Id cavity. 
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To determine the energy eigenstates of the Hamiltonian in Eq. fllSI) . \E), with the 
boundary condition in Eq. f|T9|) . one may posit 



1^) = 1 dx fE{x)^\x) + Ai,at| |0) , (20) 

where |0) is the simultaneous eigenstate of Pi and cr^ given by |(pi = — |^), i)- From 
i^Ferl^) = ^1^), one gets 

- ^^^^^ + Y\Eei'''5{x) = EMx) 

XEe{6) + YfE{0)e-i''' =E\e , (21) 
which is solved by 

fj,{x) = -^[e'ieix) + e-'ie{-x)] , (22) 



provided that 

-L + ^ + ^ = ; vt^n(^)+L-X— = . (23) 
V V 2 \v J E — e[o) 

In the inset of FigJT], Eqs. fl23p are graphically solved for e(5) = 0, using the dimensionless 
variable x = LE/v. 

For g > 1, inhomogeneities in the fabrication parameter Ej provide an irrelevant 
perturbation, since the pertinent operator scales as (^-^j ^ [9] and, thus, does not alter 
the main results of our analysis. Furthermore, today 's technology allows to fabricate 
superconducting devices with values of g ranging from g < 1, to g ^ 2 [T5]. 

Operating a YJJN near the FFP allows to engineer a realistic finite-size two-level 
quantum device with enhanced quantum coherence. Indeed, for a YJJN of finite size 
L, the FFP is stable against small fluctuations of the flux $, provided that v/L is 
sufficiently big: if 7 = 7r/3 is displaced by a small amount u, v/L needs to be larger than 
the energy splitting E^ sin(z/) between the minima of the two triangular sublattices. 
When v/L < sin(z/), there is a flow towards the SFP and, depending on sgn(z/), 
the minima of the boundary potential lie on either one of the triangular A and B 
sublattices (see FigjS]). The parameters Pi, P2 may be safely tuned to the degeneracy 
values Pi = l/(3v^), P2 = hj resorting to multipolar magnetic coils [I6] inserted in 
loops connecting the bulk superconductors at the outer boundary of the YJJN since, 
for sufficiently long chains, the magnetic flux generated by the coil does not alter the 
flux threading the circular Josephson junction array C. 

Josephson networks where n finite chains are connected to a central circular array 
C share properties similar to a YJJN. For n = 4, the resulting network is the tetrahedral 
qubit proposed in Ref . [T7] . 

In summary, our analysis of YJJNs provides an explicit example of a situation in 
which quantum impurities may be pertinently used for engineering quantum devices 
with enhanced quantum coherence. 
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Figure 4. Qualitative behavior of x"(ci;)/w in the various regimes. 
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Figure 5. Sketch of the RG flow for 7 = 7r/3 + {vj-K < 1). 
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